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We propose a simple experiment tliat everybody can carry out just while reading this paper. The 
only thing we need is a pen or a pencil and a finger. The dynamics of the falling pen, which is 
in touch with the finger, depends essentially on the initial inclination angle of the long axis of the 
pen. We point out two different types of trajectories of the falling pen and we investigate initial 
conditions for each of them to be realized. We write differential equations describing the dynamics 
of the pen and analyze them both numerically and analytically. 



I. SETUP AND EXPERIMENT 

We would like to show in our paper how sometimes 
very usual things taken from our daily live can be formu- 
lated into an elegant physical problem. Let us carry out 
a simple experiment. We take an ordinary pen by one of 
its ends while keeping the other in touch with our finger 
as it is shown schematically in Fig. ^ At first, we allow 
for a large enough initial angle between the long axis of 
the pen and the vertical, say ao = 7r/4. We let the pen 
fall down. Due to the contact with the finger, however, 
it is not a free fall. While falling down the pen slides 
across the finger and this makes it to rotate. Finally, it 
fell on one side from the finger (for large initial angles, it 
is a side where the center of mass of the pen was located 
initially). If then we decrease gradually the initial angle 
and repeat the experiment again, at a certain angle we 
notice that the center of mass of the pen starts to pass 
through the point of the "pen-finger" contact and the 
pen falls down on the other side from the finger. 

The question is: What is the critical value of the ini- 
tial inclination angle in the considered setup which does 
separate the two different types of motion of the pen in 
our experiment? 

First, let us understand upon what this critical angle, 
ttcr, could depend? After some trials we can see that it 
depends on the initial position of the contact point be- 
tween the pen and the finger. More precisely, being di- 
mensionless, the critical angle depends on the ratio of the 
initial distance between the center of mass of the pen and 
the contact point to its total length, Iq/L. Interestingly, 
this ratio together with the coefficient of kinetic friction, 
fi, between the pen and the finger, which is also dimen- 
sionless, are the only parameters which determine the 
critical angle (we will consider a pen and a finger of zero 
thickness which is close to realistic conditions and sim- 
plifies essentially the problem). Indeed, it is not difficult 
to see that we cannot construct any other dimcnsionlcss 
parameter in our problem, for example, from the mass of 
the pen or the gravitational constant, which means that 



the critical angle does not depend upon them. 

We would like to determine the critical angle for all 
values of the parameter Iq. However, this proves to be 
rather a hard task. We will solve it only numerically. 
Two limits are of particular interest in this problem that 
correspond to the case of /q = and lo ^ L/2. It is 
interesting to test our intuition and to try to guess the 
value of the critical angle in these limits. Is it 7r/2 and 
respectively or not? We will try to exploit both of these 
cases analytically. 



II. EQUATIONS OF MOTION 

Differential equations which describe the evolution of 
the solid body momentum, P, and angular momentum. 
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FIG. 1: The pen and a finger are in contact at the point O. We 
show also a coordinate system, forces, and basic parameters 
which are defined in the text. 
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M, can be written in the general formi: 



dP 

'dt 



dM 
IT 



K, 



(1) 



where F is the total external force applied to the body 
and K is the total torque. 

There are two independent variables in our problem: 
the angle a and the length I, see Fig. ^ However, the 
normal force, N, which acts on the pen at the pen-finger 
contact point and which is perpendicular to the pen long 
axis will also enter our equations of motion. So we will 
need three equations in total to describe the dynamics of 
the pen. It is convenient to derive two of these equations 
from the angular momentum evolution equation written 
with respect to two different axes. First, we write this 
equation with respect to the axis passing through the 
center of mass of the pen perpendicular to the plane of 
Fig. n (ij^ order to exclude the weight force from consid- 
eration): 



dt 



(2) 



where Iq is the proper moment of inertia of the pen, ui{t) 
its angular velocity, N{t) the absolute value of the normal 
force, l{t) the lever of this force. The dependence on time 
for quantities we indicate by the parentheses on the right 
with the time variable, i, inside. 

Now, we consider the evolution of the angular momen- 
tum defined with respect to the axis passing through the 
contact point O perpendicular to the plane of Fig. (we 
exclude the normal force from the equation in this case) : 



dt 



[{■mP{t) + Io)uj{t)] =m,gl{t) sin a{t), (3) 



where m is the mass of the pen, g the gravitational con- 
stant. 

The last equation we will obtain from the evolution 
equation for the total momentum of the pen. For this 
purpose we choose the frame of reference which rotates 
together with the pen, contact point O being in its origin. 
In this frame of reference the center of mass moves along 
one of two coordinate axes which coincides with the pen's 
long axis. Since this frame of reference is not inertial we 
should take into account the centrifugal force acting on 
the pen as well. Finally, we have: 



Then, we find from Eq. (0) the normal force: 

Nit) = AM^ 
^' l{t) dt 

and insert it in Eq. I^J which takes now the form: 

dviit) 



(6) 



dt 



= l{t)uj\t) - g cos a{t) 



mP{t)+Io 



gsma{t)-2^Lj{t) 
at 



(7) 



We note that by definition uj{t) — da/dt and vi{t) = 
dl/dt. Thus, we rewrite Eqs. Q and (jT)) and we obtain 
the system of two coupled nonlinear second order dif- 
ferential equations which describes the evolution of the 
angle and the center of mass position in time: 



I" = '2 



r' + io 

a " — cos a 



(sin a - 2l'a') 
l^ + Io 



(8) 

(sina-2/'a'). (9) 



Here we have made all quantities dimensionless upon the 
substitutions: t :— t^/g/L, I := l/L, and /q := lo/imL^). 
We use a prime instead of d/dt from now on in our for- 
mulas in order to make them more compact. 

We require that the solution of this system would 
satisfy four initial conditions: a(0) — ao, 1(0) — Iq, 
a'(0) — 0, and l'{0) ~ 0. The last two conditions corre- 
spond to the zero initial velocity of the pen. 

To find a solution of the above system is not a simple 
task. One way consists in employing numerical methods. 
In the next Section we will show results which gives this 
approach. After that we will consider two cases when the 
system can be simplified and investigated analytically. 

In Appendix^we give an alternative way of derivation 
of the above system of equations specifically for those vig- 
ilant readers who could feel themselves uncertain about 
the rotating frame of reference. 



III. NUMERICAL SIMULATION 



dt 



[mvi{t)] ^ ml{t)uj^it) - mg cos a{t) + Fi,{t), (4) 



where vi (t) is the velocity of the center of mass along the 
pen's long axis, F{^{t) — iiN{t) the friction force. 

We will now simplify the above equations. First, we 
transform Eq. © to obtain: 



duj{t) 
dt 



ml{t) 



mP{t)+Io 



gsma{t)~2'^Lo{t) 



(5) 



In order to solve numerically the system of Eqs. © 
and ((SJ we first come back to variables to and vi and we 
rewrite it as a system of four first order differential equa- 
tions with four independent variables: a, oj, I and v :— vi. 
Then, we use the finite difference method, namely Euler's 
method, to solve it. For this we divide the time interval 
[0, T] into N small intervals, each of a length At = T/N, 
and substitute, e.g., the derivative da(t)/dt by the ratio 
[a(ii+i) — a(<i)]/A<, where =ti-\-lS.t and a (to) = do- 
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FIG. 2: Trajectories of the center of mass of the pen for three 
different initial positions. 



We show in Fig. O the calculated critical angle val- 
ues for the whole range of the initial length ^o- We have 
carried out the calculation for three values of the co- 
efficient of kinetic friction /i. In the case of non zero 
value of the coefhcient we see that acr is smaller than 
7r/2 when Iq = 0. This is due to the fact that there is no 
sliding for the angles a > tt/2 — arctan/i, when the fric- 
tion is available, which follows from the couple of equa- 
tions at equilibrium: mg cos a — fiN and mgsma = N. 
Thus, when a > tt/2 — arctan/i the pen first rotates 
without sliding, its center of mass falls below the finger 
and then (when a > tt/2 + arctan/i) it starts to slide. 
In contrast, when a < tt/2 — arctan// the pen starts to 
slide from the very beginning. The numerical result for 
Iq = 0, see Fig. O agree well with the critical angle value 
Q!cr = tt/2 — arctan/i. 

A peculiar feature of the dependence of on Iq is 
its finite value for Iq — 1/2, i.e. when initially the pen 
touches the finger by its very end. This angle is about 
27° for /< = 0. 



So, starting from i — we have at each step: 



IV. ANALYTICAL SOLUTION 



At- 



k 



(sinofj - 2viUJi) , 



Vi+i ^ Vi + At 



lujf — cos ai 



k+i = li + Atvi. 



I' + Io 



(sinofj - 2viUj,) 



(10) 



We consider now two cases when the system of equa- 
tions © and lini, which describes the evolution of the 
inclination angle and the center of mass position of the 
pen in time, can be substantially simplified. 



Surprisingly, this simple numerical scheme turns out 
to be stable and works very well. We show in Fig. |3 
the results of our calculation for the trajectories of the 
pen's center of mass for three different initial positions. 
We use the value of the proper moment of inertia for a 
homogeneous stick Iq — 1/12. We note a complex be- 
havior of the pen when its center of mass passes through 
the contact point twice. In two other cases the pen falls 
explicitly on different sides from the finger. 

In order to determine the critical angle we simulate the 
fall of the pen for a given initial length ^o- First, we take 
the initial angle ao, to be equal to 7r/4 and calculate the 
trajectory of the pen. If the pen fell on the left side from 
the finger (we just look at the position of its center of 
mass at a height y = 0.1) we subtract from ao its half 
and repeat the calculation. Otherwise, if the pen fell on 
the opposite side, we add to the critical angle the half 
of it. We repeat the cycle about ten times increasing or 
decreasing the initial angle at each step by a half of an 
increment from a previous step. This procedure allows to 
achieve a very good precision of the critical angle value 
of about 1/2^" 1Q~^. The above method lies in the 
class of dichotomy methods. 




0,1 0,2 0,3 0.4 0,5 

length, Iq, in units of L 



FIG. 3: Numerical results for the critical angle for three dif- 
ferent coefficients of friction: sl) fi — 0, h) jj, — 0.5, and c) 
jj, — 1. Line d) corresponds to the analytical solution in the 
case of critical angles close to tt/2. 
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A. Angles close to 7r/2 

It should be clear that the critical angle will be about 
7r/2 when Iq is close to zero. In the limit |^| <C y/To and 
I7I <C 1, where 7 = 7r/2 — a, Eqs. ||SJ) and © simphfy to 
the system (we consider only the case = 0): 



/07" 



-I, 
-1- 



(11) 
(12) 



We have supposed the derivatives 7' and V to be negligi- 
bly small compared to unity as well. The solution of this 
system is a combination of trigonometric and hyperbolic 
sines and cosines. However, due to the zero initial con- 
dition for the velocity the sines drop out of the solution. 
Thus, we have 

^{t) — A cos ijjt + B cosh ujt^ (13) 
l{t) ^ Ay/kiCOSU)t- B^/l^coshujt, (14) 

where A and B are constants and oj = Iq . 

We suppose that the critical trajectory corresponds to 
the case when Z(T) = and 7(T) = for a particular 
moment T . This could probably contradict the observa- 
tion that the center of mass can pass twice through the 
contact point (see Fig. However, we can hope that 
our assumption is valid at least for the limit 7^1. It 
results in B = and the oscillatory solution for 7 and I: 

^{t) — ■yo cos Ljt and l{t) ~ Iq cosujt, (15) 

where 70 — Iq/^/Iq, which is exactly the condition we 
searched for. Finally, the critical angle in the considered 
limit equals to acr = ■''■/2 — Iq/^/Iq. We trace this depen- 
dence of acr upon ^0 in Fig- 121 The agreement between 
our theoretical result and numerical data for angles acr 
close to 7r/2 is quite remarkable. 



B. Small angles 

Another limit of interest corresponds to small angles a. 
The principal question is whether it is possible to prove 
analytically that for = 1/2 the critical angle acr has a 
finite non-zero value. So, we consider small angles a <C 1 
and ^0 close to 1/2 and we rewrite Eqs. © and ® in 
this limit as: 



+ Io)a" + 2ll'a' -la^Q, 



(16) 
(17) 



The solution of Eq. is I = -i^/2 + which corre- 
sponds to a free fall of the center of mass. In Eq. H16() we 
switch from the time variable, t, to the length variable, 
arriving at the equation: 



da 



21 



P 



2(1 - lo 



da 



Ic 



2(P+/o)(/o-0 



= 0. 

(18) 



We are not able to fully solve this equation. However, 
if we consider only the very beginning of the movement, 
when (Zo — Z) <C 1, we need to keep only third and fourth 
terms in Eq. (|18|) which simplifies it essentially: 



da 

m 



Ic 



= 0. 



The solution of the above equation is 



a{l) = ao\ 



P + Io' 



(19) 



(20) 



which fulfills initial conditions a = ao and a' = at 
t ^ 0. 

In order for our approximation to be self-consistent we 
should have a(l) <C 1 for all values of I. While the solu- 
tion (|20|) is valid only for (Zq — Z) <C 1 we assume it does 
not differ much from an exact solution even for I close 
to zero. We cannot prove this assumption analytically, 
though a numerical simulation of Eq. H16|) supports it. 
Thus, from the con dition of s elf-consistency for / — 
we should have aoyZf/io + 1 <C 1. This means that, 
in the case of Zq = 1/2, the angle a is always small if 
ao 1/ \/Iq/ Id -\- 1 = 1/2 (we note that this limit is 
close to 7r/6 = 30°). For these small initial angles the 
center of mass of the pen will pass across the contact 
point with the finger. In contrast, for larger initial an- 
gles the pen will fall on the side where its center of mass 
was located initially. So, the critical angle is in between 
and it is larger than zero. 



V. SUMMARY 

We have proposed a simple experiment, which needs 
a pen and a finger all in all, and we have formulated it 
into a mathematical problem. The key question of the 
problem is the value of the critical angle which separates 
two different types of trajectories of the falling pen. Nu- 
merical methods allow to solve this problem and to find 
critical angles for any initial center of mass position Iq. 
Theoretical analysis of the equations provides us the crit- 
ical angle value in the limit of small Iq and allows us to 
say that for Iq close to 1/2 the critical angle should have 
a finite non-zero value. 

We would like to refer the reader to several problems 
which are also about the dynamics of the falling pen, 
however, in a totally different setupiSii4 We have en- 
countered them while studying our problem. They are 
quite nice! 
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APPENDIX A: AN ALTERNATIVE DERIVATION 
OF EQUATIONS OF MOTION 

We keep the equation which describes the evolution of 
the angular momentum in the center of mass frame of 
reference: 



dt 



[Ioio{t)] = N{t)l{t). 



(Al) 



It permits to express the normal force through a and I 
as TV = Ioa"/l. 

We write then two equations for the evolution of the 
pen momentum in projection on x and y axes: 

mx" = N cos a — fiN sin a, (A-2) 
my" — mg — N sin a — fiN cos a. (A3) 



Next, we introduce polar coordinates according to: x = 
—/sin a and y = —I cos a. We note the expressions for 
the second order derivatives in these coordinates: 



-x" = I" sin a + 21' a' cos a + / (a" cos a — (a')^ sin a) 
- y" — I" cos a — 21' a' sin a — I (a" sin a + (a')^ cos a) 



We make variables dimensionless as we did in Sec. ^ 
and finally we obtain Eqs. © and © from Eqs. (|A2|I 
and ljA3|l rewritten in polar coordinates. 



^ L. D. Landau and E. M. Lifshitz, Theoretical Physics I: 
Mechanics, (Pergamon Press Ltd 1980). 

2 A. Chernoutsan, "Where Will a Stick Slide to?" (in Rus- 
sian), Kvant 4, 41-42, 1998 (Moscow Center for Continuous 
Mathematical Education). 

^ A. Chernoutsan, "The Stick Continues to Fall..." (in Rus- 



sian), Kvant 2, 26-27, 1999 (Moscow Center for Continuous 
Mathematical Education). 

L. Turner and J. L. Pratt, "Does a Falling Pencil Levitate?" , 
Quantum, March/ April 1998 (National Science Teachers 
Association (NSTA) and Springer- Verlag New York, Inc). 



